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We introduce a simple and straight-for-ward averaging procedure, -which is a generalization of one 
which is commonly used in electrodynamics, and sho-w that it possesses all the characteristics we 
require for linearized averaging in general relativity and cosmology - for weak-field and perturbed 
FLRW situations. In particular we demonstrate that it yields quantities which are approximately 
tensorial in these situations, and that its application to an exact FLRW metric yields another FLRW 
metric, to first-order in integrals over the local coordinates. Finally, we indicate some important 
limits of any linearized averaging procedure with respect to cosmological perturbations which are 
the result of averages over large amplitude small and intermediate scale inhomogeneities, and show 

• our averaging procedure can be approximately implemented by that of Zotov and Stoeger in these 

' cases. 

a^ ■ 

^ ' I. INTRODUCTION 

Ph. 

, It is usual in cosmology to consider the standard Universe as spatially homogeneous and isotropic on the largest 

■ scales. In fact, there is very good observational support for doing so. Ho-wever, we also kno-w that inhomogeneities 
exist at almost all scales - the smaller the scale the larger the inhomogeneity. All these inhomogeneities are apparently 
insignificant in cosmology, as long as we are not interested in modelling structure formation. It is therefore usual and 
considered acceptable by the vast majority of researchers to ignore them in investigating the dynamics and geometry 
of the Universe as a -whole. Nevertheless, Ellis 0] has given compelling reasons -why conceptually the large-scale 
cosmological metric should really be an average over very large regions of space-time. According to this vie-w, the 

■rj" ' metric g„j^ and the Einstein field equations 

^: 1 

■ G,,. = R,,. - -g^,R = 8ttT^,, (1) 

' -where R^^^ is the Ricci tensor and R is the Ricci scalar, both depending on second-order derivatives of f/^,y, and T^^ 
is the stress-energy tensor, must be averaged over small and intermediate local inhomogeneities on larger and larger 
I ' scales to obtain the average cosmological metric and the averaged dynamical equations. 

99- 

^ I Ho-wever, averaging and operating -with the Einstein differential operator on a metric do not commute, because of 
" '~j ■ the nonlinearity of the operator. Thus, the solution to the averaged Einstein equations -will not be the averaged metric. 
r> The averaged metric will, therefore, obey equations different from the averaged Einstein equations. In general, the 
I smoothing-out operation -will introduce extra tensor terms in the field equations, -which may affect the dynamics and 
the energy conditions in the averaged universe. Furthermore and just as importantly, it is a very complicated and 
unresolved issue to define an adequate averaging scheme -with the necessary properties - including the uniqueness of 
the averaged objects and their at least approximate tensorial character. This difficulty arises, because, in general, 
integrating a tensor field does not yield another tensor field in a curved space-time. 

In spite of these difficulties and uncertainties, several averaging procedures have been proposed by Isaacson 
Noonan Zotov and Stoeger and more recently by Boersma Q. Their different philosophies and results, 

together -with a comparison -with the averaging scheme we shall introduce here, are briefiy described and discussed 
below. Although far from the object of the present paper, it is also worth mentioning the ideas developed by Zalalet- 
dinov in constructing his theory of macroscopic gravity . He has used the concept of duality ~ the existence of two 
types of observers, microscopic and macroscopic - for the study of classical physical phenomena. Relying on this, 
he constructed a theory based on averaging a curved space-time itself and then determined the geometrical objects 
(metric, connection and curvature)which describe the averaged space-time This approach is non-pertubative in 
nature. 

In this investigation we define an averaging scheme, which is a straight-forward generalization of one often used 
in macroscopic electromagnetic theory carefully examine its properties, particularly for weak fields and for per- 
turbations, and show that it is also an improvement of Noonan's averaging procedure ^j^- The behavior of our 
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averaging scheme in the cases in which the hnear approximation is sufficient is acceptable. In these situations the 
noncommutabiUty of averaging and operating on the metric with the relevant differential operator is replaced by 
commutability, yielding simple, almost trivial results. However, these provide an essential reference for evaluating 
our averaging procedure and for reaching firm conclusions concerning the proper interpretation of and constraints 
on cosmological averaging in general, as well as interesting and important applications. Since much of cosmology 
involves perturbation treatments, we believe it is crucial to be clear about the averaging operator and its results in 
this context first, in order to be able to understand the cases in which nonlinear effects become important. 

Next, we discuss in detail the general properties of our averaging operator - in particular its approximately tenso- 
rial character in cosmological coordinates and for general coordinates in the weak field case, and its apparent lack of 
uniqueness with respect to local coordinates. 

Penultimately, we apply our averaging procedure to the important case of perturbed Friedmann-Lemaitre- 
Robertson- Walker (FLRW) cosmologies, showing that averaging an exact FLRW space-time yields, to first order, 
in the local coordinates, an FLRW space-time, and further that averaging a perturbed FLRW cosmology yields, as 
expected, another perturbed FLRW cosmology. We then discuss in some detail, from the averaging point of view, 
what type of inhomogeneities in an FLRW background may be considered as perturbations of an FLRW metric. Our 
linear formalism is able to handle deviations for which the local dynamics are not completely decoupled from the 
general expansion of the universe. In contrast, if the local inhomogeneities are not in the linear regime, then a first 
order perturbative approach cannot be applied. Nevertheless, our general averaging scheme is still valid, but the 
resulting metric is not a simple superposition of FLRW plus perturbations. 

Finally, we show how our averaging procedure may be approximately implemented in this nonperturbative case in 
simple situations. The procedure used by Zotov and Stoeger 1^, applied throughout a space at each point, essentially 
accomplishes in an easily implementable way what our averaging procedure requires. 



This is a simple extension of the definition used by Jackson for macroscopic electromagnetic fields. Here Q is any 
field - it may be a scalar, a vector, or a tensor. The coordinate x gives the cosmic, large-scale location of the averaging 
volume in space-time, and x' gives the small-scale location of a point within the averaging volume Vl' relative to its 
location x (its center, if it is a sphere) in the space-time, both expressed in Minkowski (rectangular) coordinates. 
Throughout this paper we often write x and x' for the cosmological and local coordinates, respectively, and x + x' 
for their sum, by which we mean x'' , x'^ and a;'' -I- a;''' , respectively. The metric of the space-time, for the observer 
who sees the inhomogeneities, is g^v{x + x') at the point x + x' . This is why the determinant of the metric g in 
the integrands of Eq. (2) depends on x -\- x' . Though each integration is only over the local regions of space-time 
dominated by the inhomogeneities, there are innumerable such averaging volumes spread over the observable universe 
- one for each cosmological coordinate x - and the metric within each one depends, of course, not just on the local 
coordinate x' but also on the cosmic location of the the averaging volume, given by the cosmological coordinate x. 
Obviously, the result of the averaging will depend strongly on the length scale over which it is performed. In general, 
that will be larger than, or of the order of, the characteristic length scale of the inhomogeneities over which we want 
to average. 

This definition is not invariant under a change of coordinates. It can only be implemented as such in a coordinate 
system like Minkowski's in which the cosmological coordinates and the local coordinates are "parallel", or transla- 
tionally related. Nevertheless, as we shall see, our averager has some very nice properties. In particular, it is almost 



II. THE AVERAGING PROCEDURE 



A. Definition 



We define an averaging operator, or simply, an averager^ acting on a field Q as 




(2) 
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invariant under transformations of the cosmological coordinates, and under transformations of the complete coordi- 
nate a; + a;' in the weak field and perturbation cases. 

We would like to stress that the quantities that are averaged are always referred to the same scale. In the process 
of averaging, the only observer involved is the one at the local 'small scale', which means that we are treating the 
problem in a self-consistent way. In other considerations one might have another observer - a cosmic observer - who 
does not see the inhomogeneities. In those cases, we would need to compare the results of the latter with the averaged 
ones. 

There are two motivations for choosing this procedure for averaging. The first is simply that it implements what 
we intuitively envision as averaging over a given length scale, providing an assignment of an average value of a quan- 
tity to every point x throughout the space-time manifold taking into consideration the possibly different values of 
the metric in different cosmological locations, and specifying unambiguously the relationship between the large-scale 
cosmological coordinate system and the local, small-scale coordinate system, over which the averages are carried out. 
One simply takes the averaging volume and shifts it from point to point throughout the universe, averaging at each 
point according to Eq. (2). This is precisely what is done in making the transition from microscopic to macroscopic 
electromagnetism [||, as we have indicated above. Secondly, as we can explicitly show (see Section IV), the averaging 
procedure defined in Eq. (2) yields quantities which are approximately tensorial in character with respect to the 
cosmological coordinates x, and with respect to the full coordinates x + x' for weak gravitational fields. That is, if 
Q{x + x') is a tensor with respect to a; -I- x' , < Q{x) > is almost a tensor of the same type with respect to x, and 
with respect to a; -f a;' for weak fields and perturbed space-times - deviating from a true tensor by only very small 
quantities. This is a very attractive property of our averager. 

At the same time, there are some other aspects of our averager which, at first sight, are not so attractive and 
require further explanation. First, as we have already mentioned, the relationship between the cosmological and the 
local coordinate systems, as formally expressed in Eq. (2), can only be realized in a very retricted class of coordinate 
systems - those in which they can be related to one another by simple translation, so that the "complete" coordinate of 
any point can be expressed simply as the sum of the cosmological coordinates and the local coordinates at that point. 
Most coordinate systems will not fulfil this requirement. This may not be a problem, since, within the demands of 
coordinate covariance, we are allowed to choose any coordinate system we want. In particular, any coordinate system 
which simplifies the formulation of our problem may be selected - as long as the quantities in question are tensors, 
or nearly tensors, thus assuring us that they are equivalent to their forms in more complicated coordinate systems. 
For other choices we could also specify the relation between the cosmological coordinates and the local coordinates, 
but the relationship between the two would, in general, be much more complicated - not specifiable by a simple sum 
of the two coordinates, that is by simple translation. 

This relation between the cosmological coordinates and the local coordinates leads another potential difficulty. 
Although it is clear that the average we have defined in Eq. (2) is approximately covariant with respect to the 
cosmological coordinates [Q, as we shall show in Section V, strictly speaking a transformation of the cosmological 
coordinates x in the integrals of Eq. (2) should be accompanied by a transformation of the x' coordinates, along 
with an induced change in their functional relationship (from a simple sum to something more complicated), as in- 
dicated above. However, even then the averaging over the local coordinates would still not be invariant under those 
changes. In light of the arguments given by Ellis and Matravers jl^ and Ellis, Matravers and Zalaletdinov for 
using preferred coordinate systems in general relativity which are appropriate to the physical situation and simplify 
the problem, even though they break coordinate invariance, we maintain that this lack of covariance of the averaging 
procedure with respect to the local coordinates should not be considered an essential problem in most cosmological 
applications. As we shall show, this is certainly the case for weak gravitational fields and for perturbations from an 
exact cosmological solution (e.g. FLRW) to the field equations. In these cases the deviations caused by the lack of 
coordinate covariance of the integrals with respect to the local coodinates are small. Even for more general cases, it is 
likely that these deviations will be small, as long as we are averaging over volumes for which the space is almost flat 
- for which the length scale of the averaging volumes is smaller than the radius of curvature of the universe. Thus, in 
that sense the most important property of our averaging procedure is that it is covariant relative to transformations 
of the cosmological coordinates. 

Another way of expressing the above objection is that in carrying out our averaging procedure we are effectively 
adding vectors and tensors at different points, which is not a well-defined, or even an allowed, operation on vectors 
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and tensors if the space-time is curved. We should really incorporate bi-vectors in the averaging integrals in order 
to translate vectors and tensors to the point designated by the cosmological metric and then add them, as proposed 
by Isaacson . The principal reason why we have not done this is that there is no easy way of defining the needed 
bi-vectors without already knowing the cosmological background metric, and it is precisely this cosmological back- 
ground metric that we eventually want to determine by carrying out the averaging over small scales. We discuss this 
in more detail in the next subsection below. (Even though throughout most of this paper we limit ourselves to cases 
in which we have a well-defined background - the weak-field and perturbed FLRW situations - we want to be able 
to use our averager in more general cases. We shall present one way of doing that approximately in the last section.) 
Furthermore, as we have just discussed, in the case of weak fields and of almost flat averaging volumes, the errors 
introduced by neglecting the bivectors will be small. Zalaletdinov ||^ constructs his macroscopic theory of gravity 
with bi-vectors, but these have to be solved consistently with the other field equations. This is a possible way of 
proceeding, but leads to such a complicated theory that solutions in the simplest cases have yet to be obtained. 



B. Comparison with other definitions 

There are a number of other definitions of averaging which have been suggested. We describe them briefly here, 
and compare them with ours, indicating its advantages. 

• Noonan's Operator 
The averaging operator given by Noonan MM is 



Q'ix') ^-g{x') dn' 

< Q' >= ^— , (3) 

-gdn' 



where 



— g refers to the determinant of a macroscopic metric. 

— The region of integration is finite and its size is bigger than the usual scale of the microscopic observer and 
smaller than that of the macroscopic observer. Using this fact the denominator may be approximated by 
^ dn in the cases in which the variation of g is not abrupt. 

— The dependence of < Q' > on the space-time coordinates comes from using indefinite integrals. This means 
that the variable in the averager is now the boundary of the integration volume - which, we believe, is not 
well-defined for an averager. 

The three points mentioned above represent the differences from our definition of the averaging procedure. One 
possible weakness of Noonan's approach (for details, see also |fj]) is the simultaneous introduction of two scales 
in his averages: the metric in his definition is one of large scale, whereas his averaging quantities are over small 
scales. 

• Isaccson's Operator 

Isaacson's procedure of averaging |^ directly deals with the problem that, in general, the result of integrating 
a tensor field does not give another tensor, because tensors at different points have transformation properties 
which depend on their location, as we have already indicated. Since one can only add tensors at the same 
point, the objective is to carry them to a certain common point and to add them there. To do that, one has to 
introduce the bi-vector of parallel displacement j^{x, x') [p^~p^. This object transforms as a vector with respect 
to coordinate transformations at x or at x' , and it has the property that, given a vector (or tensor, in general) 
Ap at x' , then Aa{x) — ja{x, x')Ai:j{x') is the unique vector at x that can be obtained by parallel displacement 
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of from x to x' along a geodesic. Q 

Given a tensor T^^ Isaacson's averaging operator is defined as: 

<T^u{x)>^ I r^{x,x')fJ{x,x')T^.0>{x')f{x,x')d^x' (4) 

J all space 

where f{x,x') is a weighting function which faUs smoothly to zero when x and x' differ by a distance greater 
than d, its integral is normalized to one over all space-time. This definition carries with it all the properties of 
the bi- vector of parallel geodesic displacement, for example 

j^' (a:, {x, x')ga'f3' {x') = g^,^{x). (5) 

This implies that if we select ga'P'ix') to be the small-scale metric, it is impossible to obtain an averaged metric 
(corresponding to cosmological scales) different from the assumed large-scale metric. That is, 

< 9fj.u{x) >= / j^^'{x,x')j^'{x,x')ga'i3'{x')f{x,x')d'^x' = 

J all space 

9iiuix) f{x,x')d'^x' = gf,^{x). (6) 

J all space 

This represents a crucial problem for averaging, since, as we have stressed, we ideally want to obtain metrics 
on larger scales as averages of the small-scale metrics, without assuming the particular form of the large-scale 
metrics. 

• Zotov and Stoeger Averaging 

Zotov and Stoeger first presented a simple procedure in which they average over an unbound distribution of stars 
-in a static background- and of galaxies -in an expanding background |^ . They use an averaging scheme which 
recalls a three dimensional version of finding a running mean with the simplest elementary volume: a sphere. 
Later, using similar techniques they construct averages over hierarchies of Swiss-cheese regions in elementary 
cosmological cells, which are the smallest volumes which are expanding with the Hubble flow. Though, they do 
not provide an averaging procedure which is specifically defined at each point of cosmological space-time, but 
rather construct averages over simple inhomogeneous configurations which they assign to entire regions, one can 
argue that over large volumes this procedure approximately gives an average at each point. In the examples 
treated, the metric (Swiss-cheese) of the inhomogeneities is considered spherically symmetric and therefore either 
Schwarzschild or FLRW. The attempt is to construct the physical large-scale background metric from averages 
over inhomogeneities which may be very large in amplitude compared to the an provisional background. These 
averages are obtained by adding the averages over each inhomogeneity to one another along with the averages 
over the provisional "background" space in between the homogeneities. Bi-vectors are not used. We shall come 
back to discuss this scheme more fully at the end of the paper, as it can be construed as an approximate 
implementation of our averaging procedure in simple non-perturbative averaging situations. 

• Buchert and Ehlers averaging of Newtonian Cosmologies 

With a similar approach to that proposed by Noonan, Buchert and Ehlers [p^ljl^ ] have focused on the averaging 
problem applied to the cosmological expansion of the Universe. They restrict themselves to the Newtonian 
approach, i.e. the analog of Friedmann's equation for the motion of self-gravitating presureless fluid. They 
propose that the spatial average of a tensor fleld A in the domain D{t) should be: 

< A >D^ i / d'^xA. (7) 

^ J D 



^This is not true if there is more than one geodesic which joins the two given points x and x' , because the transport along 
each different geodesic varies. One may avoid such problems arguing that the two points are close enough to each other as to 
permit the existence of only one geodesic joining them. 
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Since the comoving volume D{t) = a^(t), the fluid elements move on the average according to 

<e>n^^^3^. (8) 

V ao 

With this in mind, aoit) becomes the new scale factor and is shown to obey an averaged Raychaudhuri equation. 
Note that there are no spatial dependencies in the averaged quantities, only time is left as a variable. This means 
that the averaging scheme yields the same result at all positions in the synchronous gauge. 

They also propose a scheme for general relativity which includes the determinant of the spatial metric g — det gij 
in the region of interest: 

< A >D= ^ I d^x^A (9) 

^ J D 

with ^ ~ j '^^^^/9^ which resembles Noonan's definition even more closely. This average was then used by 

Russ et al. to compute the effect of inhomogeneities on the age of a flat Universe. The sources of the 
inhomogeneities were taken from the Zel'dovich approximation to second order, so that their results are valid 
in an early stage of the evolution of the flat Universe. They conclude that the effect on the age of the universe 
is very small on those scales. 

• Boersma's Averaging Procedure 

Boersma derived, from basic assumptions, a generic linearized spatial averaging operation for metric pertur- 
bations from FLRW, satisfying the condition that unperturbed FLRW is a stable fixed point of the averaging 
(that is, the averaging operation does not introduce spurious perturbations in the averaged metric). He specified 
the correspondence among points in the real spacetime, the averaged spacetime and the background spacetime, 
by the introduction of a bi-tensor density in the averaging integrals, which fulfils the same function as the 
bi- vectors in Isaacson's approach. He succeeded in deriving a general form of this bi-tensor density in terms the 
background metric and the future-directed unit vectors normal to the space-like surfaces over which the aver- 
ages are being performed. With this formalism, and using Bardcen's |19 gauge-invariant quantities Boersma 
was able to resolve the gauge problem in his averaging procedure and apply it to the constraint equations on 
spacelike hypersurfaces, which are closely related to the generalized Friedmann equation. 

III. THE WEAK FIELD LIMIT 

Using our averager defined in Eq. (2) we begin to examine its properties in a very simple, almost trivial way, by 
considering its application to gravity in the weak field limit. We do this in order to assure ourselves that our averaging 
procedure fulfils the simplest intuitive requirements and to establish some results which we shall need later in applying 
it to perturbed FLRW cosmologies. 

Let us consider a Minkowskian line clement plus small corrections. That is, 

5^12^(2^) = Vf^iyi^) + hf^^{x), \hi_,^\ < 1. (10) 

Applying our averager to this metric we have 



dn' \J~g{x + x')gf^i,{x + x') 
< <7p.(^) >- r ^ • (11) 



dQ' \/—g{x + x') 



Since in the weak field limit (see below) 



^~g{x + x') = 1 + - ?7^" + ^'), (12) 
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and considering that the product of two or more elements of the matrix of perturbations is neghgible we find, to first 
(Unear) order, 

< 9,A-) >= V,. + + + 0{h'). (13) 

Thus, in general, averaging a Minkowski space + perturbations gives another perturbed Minkowski space - a result 
which is obvious, but reassuring. It is, of course, possible in some circumstances, that the second term of Eq. (9) will 
be zero, that is, the perturbations will average out over larger scales. For instance, "overdensities" may be partially 
or completely compensated by "underdensities" . 

Furthermore, it is clear that in this weak field case - and in the case of perturbations to FLRW discussed in Section 
V - the operation of averaging will commute with constructing the Einstein field equations, since these are now linear. 
Thus, averaging will not introduce any new terms in the macroscopic field equations themselves, as it does in the 
exact case 

Let us now consider a region in space, free of gravitational sources but filled with gravitational radiation coming 
from a source situated at infinity. In these circumstances we can apply the weak field limit of Einstein's equations 

V V = 0. (14) 



Let us now suppose that we have two scales, a macroscopic and a microscopic scale, both free of sources. Then on 
both scales we will see the same behavior, i.e. on each scale we can write 

h!§, = e^e"^^'''' +C.C. (16) 

for i = 1,2, where the e/jjy are the polarization tensors in each region. The microscopic scale is associated with our 
definition and the macroscopic one ought to be compared with the averaged result. 
Recalling the definition of the average in the weak field limit, 



/ d'^x' ^—g(x + x') h^iy{x + x') 
J d'^x' ^~g(x + x') 

where g = det (7^,^ = det(?7^j/ + /ip,y). Writing Eq. (^) in matrix notation allows us to compute g using 

lndet(l + A) =trln(l + A). (18) 
by properly defining A in terms of /i^^. Thus 

tr ln(l + ^) ^ ln(l + A)„„ 

a 

= ln(l + /loo) + ln(-l + /ill) + ln(-l + /i22) + ln(-l + /133) 
- ln(l + /ioo)(-l + /^ii)(-l + /^22)(-l + h^z) 

wln(-l-/loo + /lll+/l22 + /l33), (19) 



and 



.9 = 1 + /loo - /ill - ^22 ~ /133, (20) 



-g^\ + -irh^,. (21) 



This result is independent of the choice of the gauge, and the only hypothesis used is that /i^i^ is small. Going back 
to the definition of the average, and using the property just demonstrated, we find to first order in h 
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/ d'^x' hfj,^{x + x') 



Therefore, using Eq. (p 



or 



(23) 



< V>=cie^.e^'=^" (24) 

First of all, note that this averaged quantity transforms as a tensor in the weak field limit. Second the polarization 
tensors can now be redefined as they appear multiplied by a constant which contains information about the micro- 
scopic scales of the system. Thirdly, we have found that the averaged solution (i.e. the average of the solution in one 
scale) coincides with the solution on another scale. This was expected, because we have worked in the linear theory. 
There Einstein equations are linearized, and therefore, as we have already indicated, the averaging procedure does 
not introduce extra terms in them. This means, of course, that the averages are also solutions of the Einstein weak 
field equations. 

We will now extend this analysis to see what happens in the weak field limit, when considering a region in which 
there is a source of gravitational radiation. In this case, one scale would correspond to that of the source, while the 
other could be thought of as the wave zone - that is, a scale in which the length scales are much larger than the size 
of the source. When averaging, we should compare the average metric < h^^ > with the one obtained for the wave 
zone. 

The solution of the weak field Einstein equations in the presence of a source is the metric [|o| 

gtj.u = Vp.u + h^iv (25) 

with 



Combining this Eq. with Eq. ( [T^ ) we can write (to linear order in h) 

t" - \x + r " 
\x + x" — x'\ 



< K.{x) ^ / d'x'd^x" ^''''^'^''^^J' J'^tf (27) 



If we now express the energy-momentum tensor as a Fourier integral, we can analyze each Fourier component 
separately, and integrate (or add those components) afterwards. This means that we can replace in Eq. ([2^), 
5^,(f,i) =5^,(x,w)e-*"*, i.e. 

Af^ f ^—iLj{t+t"-\x+x"-x'\) 

< V(^) >= 7^7 d^x'd^x" S^Ax',Lu) ^ + C.C. (28) 

Hereafter, we drop the complex conjugate of the quantities to make the expressions look simpler. However, they 
should be recalled at the end of the calculations, as all quantities are real. So far the only approximation made was 
to consider a weak field limit. 

The scales in the problem are three: 

1. Scale of the source, given by x' . 

2. Scale of the averaging procedure (or small scale), given by x" . 

3. Scale of the metric {large scale), given mainly by x. 
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This means that scale (1) should be comparable with scale (2), and both much smaller than scale (3). Under this 
assumption we can perform a multipole expansion, which to first order is 



{x-x') 



(29) 



Thus, we can integrate over the small scale and obtain for the spatial part 



/= Id^x": 



j\x+x" —X 



^%uj\x—x I 



\X — X 



,3 // lux 

ax e 



(30) 



47r- 



^iLj\x—x'\ ' ^ 



-RCOSLUR 



■ sin LuR 



(31) 



where R is the size of the region over which we are averaging. As we stated before, R should be comparable to the 
size of the source, and much smaller than the overall scale of the problem. Substituting in Eq. (|2^), we obtain the 
averaged metric as 



< h^^{x) > 



16nG 



\x—x I 



d-'x' 



OJ\X - 



— R COS LoR - 



— 7T sinojR 



(32) 



Applying again the multipole expansion for |a"| << \x\, finally 



< h^^{x) > = 



WttG 



1 - e 



lUJ 



— RcosujR 

LU 



— ^ sin ujR 



d-^x'e- 



(33) 



or. 



< h^^{x) e*''^ e^^(f , Lu) + C.C. 



(34) 



with 



AG AttG 1 - e- 



-1 



-RcosLoR 



■ sin wi? 



d^x'e- 



'Sf_,^{x,uj). 



(35) 



Here e^;^ is the polarization tensor that an observer in the wave zone would detect. 

Note the contribution from the integration over the time t" , (1 — e~^'^'^)/iuj, with T being the size of the time 
interval in space-time. Certainly one expects the choice of the limit of integration T to be dependent on the nature 
of the problem. Let us suppose that only one frequency uj is being emitted by a source. Therefore, a suitable choice 
for T will be ^. That is, T is basically the characteristic period of the system. Note as well that we have only taken 
into account characteristic scales of the system which are observed from the microscopic perspective. 

Notice, too, that the averaged metric is, again to first order, a tensor. And even more, it has the same form as the 
solution of Einstein's equations in the weak field limit considered in the wave zone. This was to be expected, because 
we are still dealing with the linearized equations, and thus, the averaging is a linear procedure. That is, averaging 
the equations is still equivalent to averaging the metric, just because we are working in the weak-field limit. 



IV. THE APPROXIMATE TENSORIAL CHARACTER OF THE AVERAGED QUANTITIES 

In this section we examine whether or not the quantities defined by our averaging procedure are generally, or ap- 
proximately, tensors under certain conditions. 

We first analyze the tensorial status of our averaging procedure in the completely general case. In order to have 
< > as a tensor we would need the transformation of L^^ such that 
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L't^Av + x')^ -Q^-^L^p{x + x') (36) 

under a completely general function — f^{y)- If Eq. ( |36| ) is valid, we can take the derivatives outside the integral 
of our averaging definition and find that the averaged quantities behave as tensors too. To satisfy Eq. ( ^6| ) we need 
to make a generic transformation between z = x + x' and z' — y + x' , with a fixed x' . But defining z' — h(z) in this 
way is equivalent to imposing very specific constraints on the original function /: 

z' = y + x' = f{x) +x' = h{x + x'). 

This is not possible for every conceivable /, and thus the averager is not in general a tensor. 

Even when we only transform the cosmological coordinate x, so that the averaging procedure yields almost tenso- 
rial quantities with respect to the cosmological coordinates only (the averaging integrations are only over x' , so that 
the transformations of the cosmological coordinates x can be taken through the integrals), Q there is an important 
problem. Strictly speaking, we cannot transform the cosmological coordinates without also transforming the local 
coordinates and altering the functional relationship connecting them. Furthermore, the integral over the local coor- 
dinates x' is not covariant with respect to those changes. 

Can we somehow show that the average is approximately coordinate covariant in some cases? That is indeed true 
for the weak-field and the perturbed-FLRW cases. 

Let Lf^i, be a tensor. In the weak field limit our averager applied to L^,y takes the form 

Jdn'{l + lTj-0h^f,ix + x'))L,^{x + x') 

To first order in the perturbation this becomes 

<Lf,^{x)> =^ Jdn' Lf,^{x + x') + ^ Jdn'^T]"^hafj{x + x')Lf,,.{x + x') 

- ^ Jan' L^^x + x') Jdn" 1 v^'^h^pix + x"). (38) 

where fl' — Jdfl\ is the 4- volume and a real number which only depends on scale of the averaging. 

We now employ Isaacson's averaging procedure and compare it with our own. In the Isaacson's case we have the 
averaging definition given in Eq. (4), which yields a bona fide tensor. We also know that the bi-vectors j" {x,x') 
must satisfy Eq. (5). In the weak field and perturbed space-time cases this implies that 

j^'(x,x')=<-(l/2)<(x,x'). (39) 

Thus, 

< >/= / d^x'T^,{x')f{x, x') + 0(1), (40) 



where 0(1) refers to integrals of first order in h. Thus, the first term in this equation is almost a tensor -differs from 
a tensor by small terms of order h. Then, our averaging procedure gives an integral like the first term of Eq. (|4^), 
except that it contains a y/—g instead of f{x,x') -which makes no difference in the tensorial character of the integral 
in the weak field and perturbed field cases. Thus, our procedure also yields an approximate tensor in both of these 
cases, relative to the complete coordinate x + x' . 



V. THE FLRW METRIC -|- PERTURBATIONS 

Let us now consider an FLRW space-time with perturbations and apply the averager to the perturbed metric 
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9tii^{x) = g^„FLRw{x) + h^^{x), (41) 
where the FLRW Une element in Cartesian coordinate can be written |^l[] 

~ {l + ifc[a:2 + y2 + ^2]}2- ^^^^ 

Applying the averaging, we find 



< >= , . (43) 

\J-g{x + a;') d9! 



Using Eq. (|T^) the determinant of the perturbed metric is, to first order in h^^i, 

det gy,u = det gy^uFLRW (1 + /iM'^fi'FMiy)' (44) 
and replacing this in Eq. (E^) we obtain 



{giivFLRw{x + x) + /).p^(a; + x')) \J-gFLRw{x + x') (1 + h^,^g'^^\j^^)dVl' 

< g^Ax) >= 1- • (45) 

\J-gFLRw{x + a;') (1 + h^^g'^''^j^y^)d^' 



A. What is a perturbation? 

Here we shall define more precisely what we mean by a perturbation to an FLRW metric. We will say that in 
Eq. (Ell) is a perturbation if, just as in Eq. (p^), 



I VI « 1- (46) 

Thus, if h^^ is a perturbation to FLRW in one coordinate system, then in any other coordinate system, in which we 
can always write the metric as 

g'l^uiv) ^ g'l^uFLRwiv) + K'^{y)^ (47) 

h'^j^ is also a perturbation. Though /i^^^ will naturally be different from /i^i^, it will always satisfy Eq. (^), if h^^i, 
itself does so. That is simply because, if we apply a general coordinate transformation 



to Eq. ([4l|), we trivially find that the coordinate transformation Eq. ( |4q ) also relates the h'^^ to the /i^^, and that, 
therefore, the h'^^, will be small, if the h^i, are. Thus also the determinate of the transformed metric given in Eq. ( |47| ) 
will be given by Eq. (|^, with the unprimed metric variables just replaced by the primed (transformed) metric 
variables. 

B. Averaging FLRW Perturbations 

Let us now obtain a general expression for the average of the perturbed metric. From Eq. ( |45| ) 

N 

< gtiv>=<gt,i'FLRw + h^^>= (49) 

with 
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N = J d'^x'^Z-gpLRwix + x')gf^„FLRw{x + x') + J d'^x' y^-gpLRwix + x')h^^{x + x') 

+ I d'^x' ^-gpLRwix + x')^gpLRwi^ + x')hap{x + x')g^^FLRw{x + x'), (50) 



and 



V = j d'^x' ^/-gpLRwix + x'){l + ^gpLRwi^ + x')hap{x + x')). (51) 

The denominator can be written as a product, the first term being the invariant volume for an FLRW space-time. 
Expanding the rest of the terms to first order in /i^j/ we get 



J d'^x' ^-gpLRwjx + x^gi_ivFLRw{x + x') ^ J d^x' ^ -gpLRwjx + x^h^^{x + x') 
J d'^x' y/-gFLRw{x + x') J d^x' ^J~gFLRw{x + x') 

J d'^x'y^gpLRwjx + X')yp1j^^r{x + x')hap{x + x')g y^^FLRW {x + x') 

J d'^x'^J-gFLRwix + x') 
j d'^x' \/-gFLRw{x + x')g^„FLRwix + x ) X 
jd'^x' yJ-gFLRw{x +'x^yFLR.w(^ + x')hap{x + a;') 



(52) 



We are now in a position to show that, for our averaging procedure, the average of an FLRW metric itself is, to first 
order in the local coordinates, always an FLRW metric, and further that the average metric of a perturbed FLRW 
metric is also an FLRW metric plus perturbations. These are among the features any averaging procedure should 
possess for it to be considered minimally adequate. 

First, since the length scale represented by the local coordinate x' will always be much less than that represented 
by the cosmological coordinates x, we can assume that the volume of integration, and therefore the domain of x' , will 
be small with respect to the overall space-time coordinates x. Thus, we can expand the various functions of {x + x') as 
multi-variable (because of the multiple coordinate variables, x' , y' , z' and t') Taylor series in these local coordinates, 
and truncate the series at first order. We can then write, for instance, 

/ I /^ ^ / \ , dg^^uFLRW , X ,\ /f-gx 
QnuFLRWXX + X ) K g^,uFLRW{x) ^ {x)x . (53) 



Analyzing the first term in Eq. (p2D, we find that, in the Cartesian coordinate form of the metric specified in Eq. (|4£ 
< ggQ >= 1. Here and in the discussion which follows, we designate the four terms in Eq. ( |5^ ) by the superscripts I, 
II, III, and IV, respectively. Further, using Eq. (|53|) we find also that 



^ I ^ , dgioFLRW , d'^x' ^ -gFLRw{x + x')x'^ 

OX J d'^x'y/-gFLRw(x + x') 



These results for the first term of Eq. ( |52| ) enable us to determine what the application of our averager to an exact 
FLRW metric yields. If we carry out the integration indicated in the denominator over the averaging volume given in 
the local coordinates x'^, we find that this term vanishes, since the integral is odd in a;'^. Thus, only even terms of the 
Taylor series will yield non-zero contributions to the average. Thus, to first order in the local coordinates, averaging 
an FLRW metric yields an FLRW metric. This result is trivial in the case of averaging a flat FLRW metric, but is not 
at first sight obvious for open and closed FLRW metrics. This result is somewhat important, because, as Boersma 
1^ has emphasized, we do not want the averaging procedure to introduce perturbations - we want the average of an 
FLRW metric to be an FLRW metric! 

If we now go back to Eq. (p2[), we see that, since the terms II, III and IV will all be smaller than the pure FLRW 
term I, we have the further almost trivial result that 
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< gtii- >= gtiiyFLRW + Perturbations. (55) 

Though these are expected and obvious resuhs, they are very reassuring. Our averaging procedure is physicaUy and 
mathematically consistent. If it did not fulfil these conditions, it would have to be abandoned. 



Studying terms II, 111 and IV of Eq. (|52|) in more detail, we find similarly that to first order 



/ d^x' ^-gpLRwix + x')hf_,„[x + ^ dhf_,^ .^d'^^J-gpLRwix + x')x'x 



^'""^ Jd^x'^^gp,nw{x + x') ^ V(^) + (^) Jd^^-g^,j,^i:c + x')x' ' ^^^^ 

and that 

< gii' >= - < C > (57) 

to first order in h and x' . 

Then Eq. (p3) becomes, to first order in the local coordinates. 



^ , , . , dgi^uFLRW I .Jd^X y/-gFLRwix + X')x ^ 
< 9t^Ax) > = gf.uFLRW(X) + 7^ ^^) rjA,/ , A + 

^x^ J d'^x'y/-gFLR.w{x + x') 

h '-^ ' 9h^^^ ^^ J d^x'y'-gFLRw{x + x')x'^ 
ox J d'^x'y'-gFLRwix + x') 



However, as we just saw above, the integrals in the numerators of the second and fourth terms of Eq. (]5^) vanish 
when taken over the entire averaging domain, because they are odd in the local coordinates. Thus, we obtain the 
somewhat mysterious and apparently trivial result 

< Ofj-i^ix + x') >= g^,vFLRw{x) + h^^{x), (59) 

to first order in the local coordinates, just what we began with in Eq. (pl|). It appears that the averaging, to first 
order in the local coordinates has no effect - that the in Eq. (|5^) is the same as that in Eq. (^). This requires 
some careful comment. 



Eq. (|5^) is simply the consequence of terminating the Taylor series in the local coordinates with the linear term. 
Any further precision using this approximation resides in the higher-order terms. In particular, any further correction 
to the already small h^^{x) will be smaller than it already is (that is, of higher order than either h^„(x + x') itself 
or the value of the local coordinates x'^). In general, we should expect that | < hfj_^{x + x') > \ < \hf^,y{x)\, but the 
possibility of it being < |/i^^(a;)| can only be explored by either performing the averaging without approximation, or 
taking the Taylor series expansion of the metric to higher orders in x'^ . 

It does not make any physical sense to go to higher orders to determine the distorsions averaging introduces in 
the FLRW background metric itself (in the closed and open cases - if we confine ourselves to averaging over spatial 
hypersurfaces, it is only in the closed and open cases that higher order distorsions are introduced by averaging), since 
it does not vary significantly over local- coordinate length scales; it is already smooth on all scales up to those of the 
cosmological coordinates. It is intuitively clear that averaging over an FLRW metric should give an FLRW metric. 
From one point of view, we might almost say that the FLRW part of the metric should not be averaged - it has in 
a sense already been averaged, or rather is the result of an average having been already performed, or assumed - in 
particular, an average of the density used in the field equations. 

With respect to the perturbed part of the metric hf^i,{x + x'), however, there is a good reason to carry out the 
averaging to higher precision. In practice, it may, unlike the FLRW metric itself, vary a great deal on local length 
scales. Averaging over those local variations may significantly decrease its magnitude on cosmic length scales, even 
reducing it to zero. In fact we know that initially there were perturbations on all scales above that given by photon 
diffusion (Silk damping). And at our epoch we see that there are very large density fluctuations (much larger than 
perturbations!) on all local and intermediate scales ~ density enhancements as well as voids. In both cases the per- 
turbations on cosmological scales are really the averages over the metric fluctuations on smaller scales. 
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Thus, when we write the perturbed FLRW metric as in Eq. ( p| ) without further specification, we are in a sense 
conflating two very different cases, which must be treated separately. The first is that in which the metric of the 
universe on all scales can be represented as an FLRW metric plus small deviations, which are present on a large 



range of scales. In this case, the expression given in Eq. (41) is correct on all scales. But, at the same time it can be 



restricted to a given scale of interest - for instance, the cosmological scale. In doing that, we would naturally define 
the perturbed component h^^(x) as 

h^„y{x) =< h^^{x + x') >, (60) 

where the hf^„{x-\-x') are small on all length scales and where the averaging procedure is well-defined as applied to an 
FLRW mertric plus perturbations: As long as the metric deviations from the FLRW metric on all scales are small, the 
averaging of h^y{x + x') over a volume of any size can be performed in a meaningful way, as we have indicated above. 
That is, the innumerable small metric fluctuations on local and intermediate scales can be averaged over to obtain 
those on cosmic scales. This would be the situation in the period of the universe up until density perturbations on one 
scale or another go nonlinear. It would certainly pertain to the epoch of recombination when the cosmic microwave 
background radiation (CMWBR) was last scattered, and to some hundreds of millions of years afterwards. Thus, it 
has important applications to CMWBR anisotropics. 



The second situation would be as at present, when the metric of the universe can be be represented by the perturbed 
FLRW metric in Eq. (^ij) only on cosmic length scales. On local and intermediate scales Eq. (^l|) is far from correct - 
the metric on these scales is nowhere near FLRW, deviations from FLRW being very large. Thus, the interpretation 
of Eq. (|l|) in this case is that on scales much larger than that of any nonlinear density fluctuation, the universe 
is close to FLRW, if the h^y{x) - the deviations from FLRW on that scale - are small. There will be lower limit 
on this length scale, below which this will not be true. Then this gfiuFLRw{x) + h^^{x) must be understood as the 
average of all the small-scale metrics - some of which will be very different from the cosmological "background metric 
- over a volume of cosmological length scale. On small and intermediate length scales we will not be able to define 
a background, such that the actual metric is always a perturbation (or small fluctuation) with respect to it. How is 
this g^vFLRW + h^i,{x) to be calculated from the large amplitude small and intermediate scale fluctuations? And how 
should we implement our averaging procedure in this case? This is an important question and will be treated briefly 
in the next section. 



In this second case, the averaging procedure we have outlined here is not at all adequate. But an adequate one 
must be constructed, if h^^{x) - the deviation of the large-scale metric from FLRW - is to have any meaning with 
respect to the metric fluctuations on smaller scales, many of which are not at all small. What is usually implicitly 
assumed in using h^v{x) in this context is to consider that we have averaged over the small and intermediate scale 
density inhomogeneities to obtain a large-scale density inhomogeneity which is found to be only marginally different 
from the background FLRW density and is therefore the source of a very small large scale h^y{x).. However, this 
uncritically avoids a number of important mathematical and physical issues which really have to be resolved before 
such a procedure is validated - particularly the effective field equations which are really operative in averaging over 
large amplitude inhomogeneities, including the contribution made by gravitational binding energy to the average, the 
background metric which is used in the averaging process, and the adequacy and validity of any averaging process 
itself involving large deviations from either Minkowski or FLRW ||l|,||j^Jl^ . We need to apply our averaging procedure 
in these situations and see what it yields. In some simple, idealized cases an approximate averaging procedure such 
as that suggested by by Zotov and Stoeger js) may, as we point out below, be equivalent to ours. 

From this brief discussion, we begin to appreciate how different these two cases are, from the point of view of 
averaging - and the rather different problems implicit in each of them. As we have seen, it is only the first case, 
involving perturbations on all scales which can be dealt with using the linearized procedures developed so far, or any 
procedures limited to perturbations. 

It is important to note, furthermore, that, as we have already seen, it is possible that in either of these two cases, the 
average of the perturbations or deviations over a large enough scale will yield zero, as in some of the weak-field-limit 
cases. In fact, strictly speaking, if a given FLRW background model is a genuine background, there should be a large 
enough scale over which averaging yields the exact FLRW model itself - that is, the average over the perturbations 
on such a large length scale gives zero. This means that the "positive" and "negative" perturbations need to balance 
on the largest scales, if the universe is really FLRW above a certain length scale. If this is not the case, then we must 
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choose the average FLRW phis perturbations as the real background. It is also clear, of course, that generally the 
average of the perturbations will possess much more symmetry than the original perturbations themselves. 



C. Deviations in the Hubble Parameter 

Of course, these issues of averaging will have consequences for the observational parameters we measure, for example 
the Hubble parameter H. It can be easily shown from a perturbation treatment of the field equations that, for a dust 
equation of state, the deviation of the Hubble parameter from its value in a best-fit FLRW model is just given by |2^] 

AH = -^S, (61) 

where 6 is the density contrast, that is (5 = — — —, where p is the mass-energy density and pb is the background 

Pb 

(FLRW) mass-energy density. 

Now, obviously Eq. ( |6l] ) depends on how S is calculated, and furthermore we want the time-derivative of the 
cosmological S - its average over the length scale of the universe. Averages over small or intermediate scales will 
give us only the local or intermediate deviations of the Hubble parameter from its FLRW value. Similar to our 
discussion of averages over metric fluctuations in the two general cases in the last subsection, determination of 6 and 
therefore 6 is relatively easy in the case when density fluctuations on all scales are perturbations, as long as we have 
the required data. Then we can use our procedure to do the averaging. But again, if we are faced with the second 
case, where the density inhomogeneities on small and intermediate scales are large, as at the present time, the correct 
averaging procedure of the density fluctuations is unclear. In fact, in this case, it is not certain that as a consequence 
of averaging the field equations at smaller scales over larger scales Eq. (^) holds. Not only that, but the perturbative 
treatment from which Eq. (|6l| ) is derived will not be valid in that regime. This issue will be investigated briefly in 
the next section, and more thoroughly in a subsequent paper. Furthermore, it is obvious again that the scale over 
which the averaging is done will determine the result. If we average over intermediate scales - or determine H using 
data sampling intermediate scales, scales at which the Hubble flow can not yet be recovered, instead of cosmological 
scales - then our calculation or measurement of H will be local, not cosmological. 

VI. AVERAGING IN SIMPLE NONPERTURBATIVE CASES 

When faced with averaging over large amplitude small and intermediate scale perturbations we can in principle 
still apply our averaging procedure. However, in general it is difficult to see how we could practically implement 
it directly. This is true even in simple idealized cases, where we have, for instance, strong spherically symmetric 
inhomogeneities. However, in such idealized situations, Zotov's and Stoeger's two-step procedure Q provides a useful 
and implementable approximation to ours. 

They first average over each individual spherically symmetric inhomogeneity, which is represented either as a 
Schwarzschild metric or a local FLRW metric plus a surrounding underdense annulus. This average is performed over 
the whole vacuole representing the inhomogeneity, without including any region outside it, and is therefore referred 
to the center point of the inhomogeneity. 

Now, very large regions of the universe can be ideally considered to be made up of collections of these spheri- 
cally symmetric inhomogeneities, or spherically symmetric clusters of them, separated by regions which are either 
approximately Minkowski - if the region is not expanding - or approximately FLRW, if the region is expanding. 
These function as a temporary or intermediate background metric, representing the metric far from the center of any 
inhomogeneity or cluster of inhomogeneities. Thus, the picture is very much like the Swiss-cheese model, except that 
we envision the background as temporary or auxilliary - to enable us to perform the second step in the averaging. 

The second averaging is performed over these very large regions - usually of cosmological length scale - consisting 
of many spherically symmetric inhomogeneities. It is approximated by: 
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ffM-(x) = -^[< > ViN + g^,,B{V2 - VxN)], (62) 

where < > is the average over a single spherically symnietric honiogeneity, Vi is the volume of each homogeneity, 
and N is the number of them in the full cosmological averaging volume V2. g^vB is the auxilliary background metric. 
In this above equation we have idealized all the inhomogeneities as identical in volume and mass. We can easily 
generalize to many different sizes, x is the cosmological coordinate at which the second averaging is centered. 

Thus, the second averaging can be conceived as done at each value of the cosmological coordinate x. We simply 
move the volume V2 around the universe - to all different values of x and perform the two-step averaging procedure. 
It can be easily seen that this is equivalent to what our averaging procedure involves, precisely in terms of moving 
the same volume around the universe and performing an average over it. Thus it may be considered as a way of 
approximating what our procedure would give in these situations of strong localized inhomogeneities. As Zotov and 
Stoeger point out, g^iuix) will generally not be an FLRW metric, even if it does not depend on the spatial coordinates 
- it will superficially look like FLRW, but the scale factor will have a different dependence on time than the FLRW 
metric. In general, of course, the result of the averaging will depend on the spatial cosmological coordinates. 

Only if the averaging is over a large enough cosmological volume and the universe is such that averaging over that 
volume centered at each point in a spacelike slice of the universe gives exactly the same result, will it yield a scale 
factor independent of the spatial coordinates. If the dependence of the average scale factor on the spatial coordinates 
is weak and its time dependence not too much different from FLRW, then the average metric may be represented 
by a perturbed FLRW metric g^i,uFLRw{x) + h^t^, where now h^i, is a large length-scale perturbation from a FLRW 
metric which is defined via the above averaging procedure. However, strictly speaking, we should just treat g^v{x) as 
the average metric over that length scale, since it will no longer exactly satisfy the Einstein field equations. |^ It will 
satisfy field equations which are Einstein's, with an extra term added, due to the noncommutability of averaging and 
forming the Einstein tensor from the metric. 

Is such averaging approximately tensorial? It seems that it should be, with respect to the cosmological coordinates. 
For the prescription for averaging does not depend on the coordinate system used for the cosmological coordinates 
themselves. But do significant problems arise with regard to the implied averaging over local coordinates? And is the 
Zotov-Stoeger averaging over the moving volume really approximately equivalent to what our averager would give? 
We shall investigate these questions, along with others in a subsequent paper. 



VII. CONCLUSIONS 



In this paper we have proposed an intuitively clear averaging procedure for general relativity and cosmology, which 
is an extension to that used in electromagnetism. We have shown that it gives approximately unique and tensorial 
results in weak field and perturbed FLRW cases, and does not lead to any significant unacceptable results in these 
cases. Furthermore, it promises to be easily applicable in cases where fluctuations on all scales are perturbations, 
such as up to the epoch in which density perturbations begin to go nonlinear. 

Finally, we explore the limits of averaging perturbed FLRW universes, indicating that averaging over very large 
small and intermediate scale inhomogeneities, in order to recover the average perturbation on large scales, requires 
applying our averaging procedure beyond the simple FLRW plus perturbations case. We show how that can be 
done approximately in simple idealized cases involving spherically symmetric objects and clusters of objects in an 
intermediate background using the Zotov and Stoeger approach. 
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